Semiclassical Klein Tunneling with Berry Phase Corrections in Graphene 
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We study the transport dynamics of semiclassical electrons in graphene in the Klein tunneling 
regime by scattering wave packets off a potential step while taking the full graphene Hamiltonian 
into account. Besides establishing the basic transmission characteristics for zigzag and armchair step 
edges, our numerical simulation suggests that the wave packet dynamics is substantially affected by 
the Berry curvature, which induces lateral shifts near the step edge reminiscent of the Hall effect of 
light. This anomalous correction can be relevant for recent Klein tunneling experiments. 

PACS numbers: 72.80.Vp 



Introduction — Graphene possesses extraordinary elec- 
tronic transport properties due to the two special points 
in the band structure where conduction and valence band 
meet. At these so-called Dirac points the energy disper- 
sion is conical, and electrons circling around each Dirac 
point acquire a nontrivial Berry phase of ±7r [![. Such a 
nontrivial Berry phase can lead to remarkable transport 
phenomena as displayed by the half-integer quantum Hall 
effect in graphene [2j-|4| and the intrinsic anomalous Hall 
effect in ferromagnetic materials [f| . The Klein paradox, 
too, is closely associated with the conical Dirac disper- 
sion and may just as well be affected by a nontrivial band 
topology. The Klein paradox refers to the almost cer- 
tain transmission of relativistic Dirac fermions across a 
potential barrier regardless of its height and width [f|. 
Considerable effort has recently been put into observing 
and exploring the Klein phenomenon in gra phe ne devices 
with evidence recently reported in [if] ] based on 
resistance measurements across a potential barrier and 
comparison with theoretical predictions for ballistic and 
diffusive transport. The authors, however, also pointed 
out unresolved discrepancies with the expected angle de- 
pendence of Klein tunneling. 

Theoretical treatments of Klein tunneling in graphene 
usually rely on single valley linearized Dirac bands and 
plane wave approaches 14 -I?}- A semiclassical wave 
packet approach, however, can give further insight into 
the effective dynamics of Bloch electrons in bands with 
nontrivial topologies [3, Ojj]. This is evident from the 
semiclassical equations of motion for Bloch electrons 

r = V k f k -kxfi k , k=-e(E + rxB), (1) 

where r stands for the electron position, for the band 
energy at crystal momentum k, e for the elementary 
charge, E and B for electric and magnetic field, respec- 
tively, while, crucially, the first cross product term rep- 
resents a correction due to the Berry curvature fik and a 
driving force k. Such a modification applies even outside 
of electron systems. For example, a Hall effect of light 



was proposed, where an electromagnetic wave packet 
scatters at an interface between two optically different 
media [2^, 21 1. While geometric ray optics may prop- 
erly describe the macroscopic classical behavior, a proper 
derivation and analysis of the electromagnetic equations 
of motion using optical wave packets revealed Berry cur- 
vature corrections that laterally shift the reflected and 
the transmitted light rays at the interface. 

In the present paper we study the transport dynamics 
of semiclassical electrons in graphene in the Klein tun- 
neling regime by scattering wave packets off a potential 
step while taking the full graphene Hamiltonian into ac- 
count. Our main findings are two- fold. First, our numer- 
ical simulation reveals in some cases unusual side shifts 
and/or kinks in the wave packet motion near the step 
edge. These features indicate deviations from a simple 
reflection law or S noil's law and are reminiscent of the 
optical Hall effect, suggesting the presence of a Berry 
curvature velocity where the driving force is the potential 
gradient at the step edge [see Eq. ([TJ]. This can be seen 
in Fig. Q] (a), which shows the trajectories of three wave 
packets propagating towards a potential step. Each wave 
packet splits into a transmitted and a reflected part upon 
striking the potential step, which is much taller than the 
wave packets' energies (i.e., the setup represents a pn 
junction). Looking at the red pair of trajectories (full 
circles), the path of the wave packet fraction in region 
I, which denotes the lower part of the step at zero po- 
tential height, clearly exhibits a shift roughly parallel to 
the step edge. A shift or kink in the wave packet motion 
can also occur on the transmission side, which is marked 
as region II and forms the upper part of the step. Sec- 
ond, we establish the basic transmission characteristics 
that, in contrast to usual Klein paradox studies involving 
graphene, include intervalley coupling, finite bandwidth 
effects and band anisotropy effects. Our results in par- 
ticular reveal striking differences in the Klein tunneling 
behavior at zigzag and armchair steps. All in all, the 
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(a) Small wave packets, k Q near K' (b) Large wave packets, k Q near K' (c) Small wave packets, k Q near F 
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FIG. 1. (Color online) Trajectories of wave packets scattering at a zigzag potential step roughly twice taller than their incoming 
kinetic energies. The zigzag step edge is marked by the solid brown line, separating the lower part of the step (region I, white) 
from the upper part (region I, grey), while the dashed black line lies normal to the step edge and serves as a guide to the eye. 
The size of the symbols along the trajectories stand for the probability weights Pi/n(r), while the black line at the bottom 
indicates the initial wave packet diameter (2oy). In (a) small wave packets with ko near K' break up into a reflected (region 
I) and a transmitted fraction (region II) upon approaching the step edge (the green arrows representatively illustrate the 
scattering process). Remarkably, the set of trajectories appears completely asymmetric with respect to the step normal, while 
the red (full circles) and the blue (full triangles) trajectories feature unusual lateral shifts and/or kinks close to the step edge. 
These features gradually disappear (b) with increasing wave packet size or (c) by moving ko closer to the T point. Note that 
the bending of the green trajectory (full squares) near the x-axis in panel (b) is a finite size effect. Please see main text and 
[U for details. 



dynamics near the step edge sensitively depend on the 
wave packet and the step potential parameters and will 
be discussed in more detail below. 

Semiclas steal wave packet propagation — To study the 
single electron dynamics we consider a spinless electron 
wave packet moving in a two-dimensional honeycomb lat- 
tice towards a sharp potential step. Although it is exper- 
imentally challenging to fabricate graphene devices with 
thin junction interfaces, a sharp potential step is the sim- 
plest setup and an instructive limit to investigate semi- 
classical Klein tunneling. 

The time evolution of the wave packet \^{t)) is 
governed by the time-dependent Schroedinger equation 
id T \ip(r)) = (H + H step )\ijj(T)), which contains the usual 
honeycomb tight-binding Hamiltonian Ho with nearest- 
neighbor hopping t [l[ and a step potential 



H, 



step 



ala r 



V(r + S)blb r 



(2) 



characterized by step height Vq (i.e. V(r) = (Vq) for 
r G I (II)). Here aj (6J) creates an electron on sublattice 
A (B) in unit cell r, while 8 denotes the intra-unit cell 
nearest-neighbor vector. Below we express all distances 
in units of the nearest-neighbor distance a, momenta in 
units of 1/a and energies in units of t. 

The wave packet evolution then follows straight- 
forwardly from a fourth order split-operator 
(Suzuki- Trotter) approximation scheme IVK 1 ")) ~ 

t/N and N is the 
For simplicity, we 



TLn=ie \iP( T = 0)), where At 

total number of iterations [2 



23|. 



impose periodic boundary conditions and use sufficiently 
large system sizes and small enough time steps to reduce 
finite size effects. 

The wave packet evolution in graphene depends sen- 
sitively on the initial wave packet setup. Although we 
evolve the wave packet numerically using the real space 
representation, i.e. 



Mr)) = Et^ r )4 + *s(r;r)6t]|0), 



(3) 



where ^&A(B)(. r \ t) stands for the wave function on sub- 
lattice A (B), it is convenient to specify the initial wave 
packet at r = using the eigenbasis of the Hq. Thus 
in the following, if not stated otherwise, the initial wave 
packet is given as a Gaussian-like superposition of posi- 
tive energy eigenstates centered close to the K' point 



Wr = 0))=E*ko,a k) ro(kKi|0), 



(4) 



where creates a state with momentum k in the elec- 
tron band Q. The momentum distribution \&k ,<Tfc,r s ) 
is centered at ko and characterized by a width cr^, while 
the position in real space is given by ro. Initializing the 
wave packet using states from both electron and hole 
bands would complicate the propagation dynamics as ad- 
ditional interference phenomena such as trembling mo- 



tion (Zitterbewegung) may appear [2 

Possible Berry curvature effect — Let us investigate 
the wave packet dynamics in graphene in more detail. 
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Panel (a) in Fig. Q] shows three sets of wave packet tra- 
jectories where the incoming wave packets emanate from 
the same point in region I, but have different central mo- 
menta k located close to the K' point (for precise mo- 
mentum values see [29[ ) . These wave packets have similar 
total energies (eo ~ OAt > 0) and sizes (set by the spread 
in momentum space <7k = 0.2/a, which corresponds to a 
real space standard deviation of ay ps 13.48a), but strike 
the zigzag step edge at different angles. At the step edge 
(Vo = 0.8t) the incoming wave packets are partially re- 
flected and partially transmitted. To follow the reflection 
and transmission dynamics, the trajectories in region I 
and II arc defined via the centre of mass motion of the 
wave packet fractions. The total probability weight of 
each fraction is given by 



P 1 (t)=J2\*a(v ) t) 



r+5ie/ 



* B (r;r)| 2 (5) 



and -Pii(t) = 1 — Pj(t), respectively, and indicated by the 
size of the full symbols accompanying the trajectories. 

To extract the semiclassical dynamics, panel (b) shows 
the evolution of large wave packets (oy. = 0.02/a, ay w 
134.8a) with similar momenta and energies as in (a), 
while panel (c) displays the propagation of small wave 
packets (ay = 0.2/a) as in (a), which, however, are ini- 
tially composed of hole states with momenta near the F 
point instead of the K' point. Here, the wave packets' 
initial kinetic energy and the step height measured from 
the bottom of the hole band are eo ~ 0.2t and Vb = OAt, 
respectively. All wave packet and step setups in Fig. 
Q] correspond to a pn junction (0 < eo < Vq). Please 
see also the supplemental material for animations corre- 
sponding to the red (full circles) wave packet trajectories 
in (a)-(c). 

There are several interesting features associated with 
the wave packet motion. Obviously, reflection and trans- 
mission of small wave packets do not simply follow a 
reflection law (9- ln = re a) or a variation of Snell's law 
(rii n sin0i n = ntransSin#trans)j respectively As noted be- 
fore, the wave packet portion in region I near the step 
edge experiences in some cases a side shift roughly par- 
allel to the edge and/or abrupt changes in the propa- 
gation direction. On the transmission side (region II), 
the trajectories can also exhibit unusual features (side 
shifts and kinks) close to the step edge, depending on 
the incoming wave packet direction and the step height, 
and can vary strongly in orientation and the transmitted 
probability weight. This is shown in the supplemental 
material, where similar simulation results (including lat- 
eral shifts and offsets) for other step configurations such 
as pp junctions (0 < Vq < eo) and armchair step edges 
are presented. Considering all trajectories together, the 
scattering process in Fig. [1] moreover appears strikingly 
asymmetric with respect to the step normal, which is 
characteristic of the zigzag edge and in stark contrast to 
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FIG. 2. (Color online) Side shifts [panels (a)-(c)] for the three 
different step scenarios presented in Fig. [T] Panel (d) illus- 
trates the definition used here for the lateral shift As and the 
incoming angle 0j n . A substantial lateral shift occurs only for 
small wave packets with initial momenta near K', supporting 
the notion of a Berry phase effect. Please see main text for 
details. 



the highly symmetric trajectories appearing at an arm- 
chair edge. 

One can attempt to quantify the overall lateral shifts 
by measuring the difference between the intersections of 
the asymptotes of the incoming and the reflected wave 
packet paths in region I with the step edge. The separa- 
tion of the intersections As = Si — So| as a function of the 
incoming angle 4>- ln is plotted in Fig. [2j This way of char- 
acterizing the side shift works best for wave packet paths 
close to the step normal where changes in the propaga- 
tion direction are rather abrupt and unaffected by finite 
size effects, which tend to occur for larger incoming an- 
gles and larger wave packets. 

Figures [2] (a), (b) indicate a drastic difference between 
small and large wave packet dynamics. With increasing 
wave packet size the side shifts gradually disappear. By 
analogy to the Hall effect of light, the trajectory anoma- 
lies therefore may arise from a Berry curvature velocity 
[see Eq. ([T])]. which is known to modify the semiclassical 
equations of motion of Bloch electrons in systems with 
broken inversion or time-reversal symmetry [l8l [lij . A 
Berry curvature term can only be effective for a finite dis- 
tribution of wave vectors and hence may be diminished or 
absent in the limit of very large wave packets and plane 
wave states. 

Since Berry curvature effects are usually strongest near 
the Dirac points, one can further check the dynamics of 
small wave packets that are initially composed of hole 
states near the T point as shown in Fig. Q] (c). Such 
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FIG. 3. (Color online) Transmission profiles for zigzag and 
armchair step edges with various step heights Vb as a function 
of incoming angle <f>in [see also Fig. [2](d)]. All wave packet 
energies and sizes are fixed at eo ~ 0.4t and Ok — 0.2/ a 
(small wave packets), respectively, as in Fig. Q](a). We note 
that with increasing step height, the finite bandwidth of the 
graphene dispersion suppresses the transmission probability 
T (see black curves). 



wave packets resemble quantum particles in an ordinary 
parabolic band and undergo a total reflection at tall po- 
tential steps and are partially transmitted and reflected 
at sufficiently low potential steps. Transverse shifts or 
kinks in the trajectories are completely absent [see also 
Fig. [2] (c)], which is consistent with a vanishing Berry 
curvature near the T point. Given these observations, 
one may therefore speculate that the magnitude of the 
side shift is roughly determined by the Berry curvature 
flux through the footprint of the wave packet distribution 
in momentum space (or by the difference thereof between 
the initial and the final states). 

Finally, we briefly comment on the role of edges states, 
which have a topological origin 3(| 31 [ . Edge states play 
only a minor role since the trajectory anomalies are as- 
sociated with both the zigzag and the armchair edge (see 
supplemental material). Furthermore, the zigzag edge 
states are located energetically just above (> Vb) and 
just below the step (< 0) and thus are energetically far 
removed for the cases considered here. It may be inter- 
esting though to consider situations where the zero mode 
edge states are relevant. 

Basic transmission characteristics — Let us now fo- 
cus on the basic Klein tunneling transmission charac- 
teristics in graphene. The easiest way to distinguish a 
zigzag from an armchair step is by analyzing the trans- 
mission profiles. Approximating the transmission prob- 
ability T = Pu(t — > oo) by the probability weight in 
region II at the largest simulation time available, differ- 
ences in the scattering dynamics are clearly visible in 
Fig. [3] Although the transmission curves in both cases 
are strongly angle dependent, the armchair curves are 
highly symmetric about <fio = 0, while the zigzag coun- 
terparts are not. Another striking feature revealed in Fig. 
[3]is that the transmission probability is substantially sup- 
pressed for tall step heights (eq <C Vb), see e.g., the solid 



black (Vb = 2.0t) and the dashed blue curves (Vb = 1.2t). 
This is in contrast to the usual Klein paradox of relativis- 
tic Dirac fermions, where transmission over a wide angle 
range approaches unity with increasing step height. The 
non-monotonous behavior of the transmission curves at 
Vb = eo, however, remains preserved (see dashed blue, 
dash-dotted red, and dash-double-dotted green curves) 
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Both the symmetry difference between zigzag and arm- 
chair step edges and the decreased wave packet transmis- 
sion for tall steps result from taking naturally the full 
graphene dispersion into account in the numerical simu- 
lation. The graphene band structure exhibits a trigonal 
symmetry away from the Dirac points and possesses a 
finite bandwidth. The latter leads to the reduced trans- 
mission probability for large potential steps due to the 
lack of appropriate scattering states on the high step side. 
The band anisotropy, on the other hand, directly affects 
the reflection and transmission trajectories further away 
from the step edge and hence is encoded in the trans- 
mission profiles. For example, a zigzag edge involves 
only intra- valley scattering. The dominant reflection and 
transmission states in the present study therefore are lim- 
ited via energy conservation and momentum conserva- 
tion parallel to the edge direction to states around the 
K' point, which, depending on the step height, can even 
give rise to a preferred overall transmission direction (not 
shown). In contrast, an armchair step potential connects 
the two inequivalcnt Dirac points, which arc related by 
inversion and time-reversal symmetry and thus produces 
more symmetric wave packet motion as a function of in- 
coming angle. 

Conclusion — We have investigated the Klein tunnel- 
ing of electrons at zigzag and armchair step potentials 
employing the full graphene Hamiltonian and a semiclas- 
sical wave packet approach. Apart from establishing dis- 
tinct transmission characteristics for zigzag and armchair 
steps, we have identified unusual wave packet dynam- 
ics near the step edge. These unusual features appear 
for both zigzag and armchair step orientations and sup- 
port the notion of Berry curvature influenced semiclas- 
sical Klein tunneling. These anomalies may help resolve 
the discrepancy found in the angle dependence of Klein 
tunneling reported in [loj . Further extension of the nu- 
merical simulations may include a magnetic field, smooth 
potential steps, which may lead to mixing with Bloch os- 
cillation phenomena, and potential barriers, which are 
useful for electron focussing and lensing in graphene [32| . 
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FIG. 4. (Color online) Wave packet trajectories complemen- 
tary to those in Fig. Q]of the main text. The wave packets' 
incoming kinetic energies are roughly eo m OAt and their sizes 
are given by = 0.2/a (ay w 13.48a) as in Fig. Q](a). in the 
two upper panels, the step orientation is parallel to the zigzag 
direction. Panel (a) contains more trajectories, but is other- 
wise the same as Fig. Q] (a), where the step height amounts 
to Vo = 0.8* (pn junction). In panel (b) the potential jump is 
smaller than the wave packet energy Vo = 0.2t < eo (nn junc- 
tion), representing the classical transmission regime. Panels 
(c) and (d) are the armchair counterparts of (a) and (b), re- 
spectively. 



SUPPLEMENTAL MATERIAL 



Wave packet animations — A set of movies corre- 
sponding to the red (full circles) wave packet trajectories 
in Fig. [T]is available at youtube. 



Armchair steps and more — Figure 2] shows a selec- 
tion of wave packet trajectories for different step config- 
urations. The step setup in panel (a) is identical to that 
of Fig. Q] (a), but here were present a wider selection of 
wave packet trajectories with varying incoming angles. 
The overall shift of the reflected trajectories is clearly 
noticeable. In contrast to (a), panel (b) shows the tra- 
jectories of the same initial wave packets scattering at a 
lower step (eo > Vo). In this case, there are not only 
lateral shifts of the reflected wave packet fractions, but 
also overall lateral offsets of the transmission trajecto- 
ries with respect to the incoming direction for larger </>j n , 
see e.g. the cyan colored trajectories (leftmost trajectory 
pair) in (b). Unusual features in the wave packet mo- 
tion also occur at armchair steps as displayed in panels 
(c) and (d) for high and low step potentials, respectively. 
The armchair case therefore produces not only trajecto- 
ries that are highly symmetric about the step normal, 
but also unusual due to overall lateral shifts and offsets. 



